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Abstract 

To locate noise sources in high-speed jets, the sound 
pressure fluctuations p', measured at far field locations, 
were correlated with each of density p, axial velocity u, 
radial velocity v, puu and pvv fluctuations measured 
from various points in jet plumes. The experiments 
followed the cause-and-effect method of sound source 
identification, where the cross-correlation coefficients, 

<p; p >, <puu; p > etc., could be related to various source 
ter ms of Lighthill’s equation. Detailed correlation 
surveys were conducted in three fully expanded, unheated 
plumes of Mach number 0.95, 1.4 and 1.8. The velocity 
and density fluctuations were measured simultaneously 
using a recently developed, non-intrusive, point 
measurement technique based on molecular Rayleigh 
scattering (Seasholtz, Panda & Elam, AIAA paper no 
2002-0827). The technique uses a continuous wave, 
narrow line-width laser, Fabry-Perot interferometer and 
photon counting electronics. Light scattered by air 
molecules from a point on the laser beam was collected 
and spectrally resolved by a Fabry-Perot Interferometer. 
To determine Doppler shift caused by air flow, the image, 
formed after the interferometer, was split into two 
concentric parts and the intensity-ratio was measured by a 
pair of photomultiplier tubes. The change in the intensity 
ratio from that created by incident laser light provided a 
measure of a velocity component. Photo-electron 
counting over short-duration, contiguous bins provided a 
time history of velocity variation u(t), v(t). In addition, a 
part of the Rayleigh scattered light was measured 
directly, without passing through the interferometer, 
using a third photomultiplier tube to obtain a time history 
of density fluctuations p(t); and finally, multiplications of 
the time series data provided puu(t) and pvv(t). Two 
separate collection arrangements were used to measure 
Doppler shifts from u and v velocity components. Fourier 
transforms of the time series data provided respective 
spectra. It was observed that the density spectra S p were 
in general similar to the axial velocity spectra while the 
radial velocity spectra S v were somewhat different. The 
p-u cross-spectra show progressively decreasing 
correlation with increasing frequency. To determine 
sources of sound pressure fluctuations microphone 


signals from 50 nozzle diameters and at polar angles from 
30° to 90° to the jet axis, were cross-correlated with 
individual flow variables. The sound pressure 
fluctuations at 30° to the jet axis provided the highest 
correlation coefficients with flow fluctuations. With an 
increase in microphone polar angle, the correlation 
coefficients decreased sharply, and beyond about 60° all 
correlation mostly fell below the experimental noise 
floor. Among all turbulent fluctuations <puu; p > 
correlations showed the highest values. Interestingly, <p; 
p>, in all respects, were very similar to <puu; p >. Both 
radial velocity v and pvv fluctuations correlated very 
poorly, except for the 30° microphone locations. By 
moving the laser probe at various locations in the jet it 
was found that the strongest noise source lies downstream 
of the end of the potential core and extends many 
diameters beyond. Correlation measurement from 
turbulent fluctuations along the lip shear layer showed a 
Mach number dependency. While significant correlations 
were measured in Mach 1.8 jet, values were mostly 
below the noise floor for subsonic Mach 0.95 jet. Various 
additional analyses, such as cross-spectra calculations, 
showed that fluctuations from large organized structures 
mostly contributed to the measured correlation, while that 
from small scale structures fell below the noise floor. 


I. INTRODUCTION 


The last six decades of research has not produced a 
unanimously accepted answer to the simple question of 
what produces noise from a jet flow. A vast number of 
earlier and current researchers have relied upon the 
acoustic analogy framework of Lighthill (1954), Lilley 
(1972) and others. Lighthill’s equation is: 

a 2 P 3 2 t, 


dr 


2 \-7 2 

— a 0 V p 


dx, dx ; ’ 


( 1 ) 


Ty = P Uj Uj + 5 y (p - a 0 2 p) 


where p is air density, p pressure, do ambient speed of 
sound, Uj, Uj velocity vectors and Ty are the elements of 
the stress tensor. In the last few years there has been a 
steady growth of opposition to such an answer, based on 
multiple issues (Fedorchenko 2000, Tam 2001). A critical 
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factor that has fueled the controversy is a lack of 
experimental verification of Lighthill sources, besides the 
U 8 dependence of sound intensity, which only confirms 
proper dimensional scaling. Experimentally, thus far it 
has been impossible to directly measure the stress 
tensors. Even if one is able to measure all possible 
turbulence statistics, the next step: which parts of the 
turbulent fluctuations actually radiate in the far field, 
becomes even more difficult to determine. The common 
practice of using microphone arrays, elliptical and 
spherical mirrors in the far field to determine source 
location have both advantages and serious drawbacks. 

The advantage lies in the simplicity and the avoidance of 
measuring the complex turbulent flow, which is modeled 
simply as a distribution of monopoles. However, the 
noise sources are not monopoles and the presumed 
distribution functions required for data analysis have 
never been verified. Sound waves undergo a large 
amount of scattering before arriving at the far field; a 
linear extrapolation of sound path may lead to a wrong 
conclusion. There is a need to determine noise sources by 
independent means, which is the motivation for the 
present work. It can be argued that such an independent 
means can neither be achieved by only turbulence 
measurement/ simulation, nor by sole observations from 
far field, but through a simultaneous measurement of 
turbulent fluctuations (cause) and the far field noise 
(effect). 

The direct correlation between the cause and effect 
was originally proposed by Siddon (see Siddon 1973, Rackl 
1973) to locate sound sources, and during the 1970’s many 
experimentalists had taken on the method. Usually, velocity 
(Lee & Ribner 1972, Seiner & Reethof 1974, Schaffer 
1979, Richarz 1980 and others) or pressure fluctuations 
(Hurdle, Meecham & Hodder 1974, Armstrong et al 1977 
and others) in the jet were correlated with the sound 
pressure fluctuation measured by a fixed microphone. One 
attractive feature of this method is that the effects of 
scattering, absoiption and refraction on sound radiation are 
automatically included by virtue of simultaneously 
extracting information from both the flow and acoustic 
fields. Another advantage in calculating the Lighthill source 
strength is that the fourth order derivative of the stress 
tensor is reduced to a second order one. Unfortunately, the 
promises were mostly unfulfilled in these earlier works, 
primarily due to a lack of non-intrusive turbulence 
measurement technique. The noise produced by intrusive 
hot-wire probes (Seiner 1974, and others) or microphones 
(Hurdle, Meecham & Hodder, 1974 and others) placed 
inside the flow contributed to most of the correlation. Later 
on Laser Doppler Velocimetry (LDV) was used by 
Schaffer (1979) and Richarz (1980) among others. Schaffer 
wrote down the vast number of approximations and 
assumptions needed to relate the experimental data to 
theory. Nevertheless, the important issue of accuracy in 


velocity spectra measurements using LDV has lingered. 

The issues of seed particles following turbulent eddies, 
various biasing errors, and a limited dynamic range makes 
LDV a difficult tool to measure turbulence spectrum. An 
approximation, used by earlier researchers in evaluating 
Lighthill’ s stress tensors, is that the contribution from the 
density fluctuation is negligible: pU.U, ~ pU.Uj, where p 
is time averaged density. Our recent work (Panda & 
Seasholtz, 2002), on the other hand, shows that density 
fluctuations have significant correlation with far-field noise. 
It is to be noted that the role of density in Lighthill’s 
formulation is ambivalent, every term in cause (right hand 
side of equation 1) and in effect (left hand side) contains a 
dependence on density. 

Recently, a molecular Rayleigh scattering based 
technique has been advanced to simultaneously measure 
density and velocity fluctuations spectra in high-speed 
flow for the first time (Seasholtz, Panda & Elam, 2001, 
2002). Fluctuations occurring over frequency range of 0 to 
50 kHz have been measured. The technique is based on 
laser light scattering from the gas molecules present in air. 
Since neither any seed particles nor any intrusive probes are 
used, the technique is free from various problems faced by 
previous researchers. Air density fluctuation is the easiest to 
measure and has the highest accuracy. An in depth study of 
correlation between flow density fluctuations and sound 
pressure fluctuations from the peak noise emission 
direction had been reported earlier (Panda & Seasholtz, 
2002). Correlation study with the radial v velocity 
component and pvv fluctuations had been reported by 
Panda, Seasholtz & Elam (2000a). The present work is a 
continuation and banks on the additional capability of 
measuring axial component (u) of velocity fluctuations. 


Causality relation: Following is a recollection of the 
principles behind the Causality relation. The free space 
solution for Lighthill’s wave equation for a field point Xf 
from distributed source points X s (figure 1) is: 


p'(X f ;t) 



’( 2 ) 


Where, p represents pressure fluctuations and the terms 
inside the integral are calculated at a retarded time to 
account for propagation from the source point to the field 
point: t 1 = t - To, To = |X r X s |/ao. The integral is taken over 
the whole jet volume V. Equation 1 shows that the stress 


a 2 (pu i u i ) 

tensor T;, has two terms. The first term has 

aXjSXj 


9 different components, out of which the contribution 
from pvv (v: radial component of velocity) and puu (u: 
axial component) fluctuations are measured in the present 
work. In the following, application of the causality 
relationship to the puu term is outlined first. Various 
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issues involved in the interpretation are discussed next. 
Finally an extension of the causality relation to the 
second term of the stress tensor, (p - ao 2 p), is outlined. 
Following Proudman’s analysis (1952), the double space 
divergence can be converted into double time derivative 
under two restrictions: (a) a scalar component of the 
stress tensor is measured along the direction of 
observation from the source, and (b) the field point is far 
enough to be in the radiation field of all sources. Under 
these conditions, the acoustic pressure at the field point 
can be written as: 

p(X f ;t) = 1 f^[puu(X s ;t / )]d 3 X s (3) 

4 7i a 0 r ^dt 

where r is radial distance. The scalar components of the 
stress tensor represent longitudinal quadrupoles, made by 
on-axis positioning of two opposite dipoles. Such a 
quadrupole has very strong directivity and, in the absence 
of convective amplification, the acoustic intensity falls as 
the cos 4 from the peak radiation direction. The peak 
radiation directions associated with 

d 2 (puu) d 2 (pvv) . 

quadrupoles, — - — — and — - — — , are respectively 

3t a^d t 

along the jet axis and 90° to the axis. The former is 
affected by refraction and convective amplification while 
the latter is not. 


T o obtain acoustic intensity an autocorrelation function 
for the above equation needs to be worked out. Usually, 
the auto-correlation function is created by multiplying 
space-time separated stress tensors separated in space by 
^ and in time by X. 



In the “causality” relationship, however; this is 
accomplished by multiplying the source integral with the 
far field sound pressure: 


(p 2 )(X f ;x) = 


( - r ^f^[puu(X s ;t / +x)]d 3 X s [p'(X f ;t)] 

\|_47t r J v dt J t 

R puu;p' (X, , x f ; x) = (p uu (X s ;t + T - x 0 ) p'(X f ; t)> 


The < > bracketed expressions imply time averages. Also, 
the correlation function R puu;p - has to be calculated after 
shifting the puu data by the propagation time x=-x 0 , or 
inversely the pressure fluctuation data by x = Xq. In 


essence, equation 5 expresses the radiated acoustic field 
in terms of a time delayed integral taken over the entire 
source volume, and consists of correlation functions 
between far field sound pressure p’ and source field puu. 

It is interesting to point out a particular advantage 
of equation 5 over 4. The Lighthill’s equation by itself is 
unable to separate hydrodynamic and acoustic 
fluctuations. The former is the pressure field associated 
with turbulent fluctuations, and the latter is the radiated 
fluctuations that propagate far away from the turbulent 
motion. Lighthill’s equation is nothing but a 
reformulation of the mass and momentum conservation 
equations and therefore every solution upholds these two 
physical laws. It is known that only a small part of the 
disturbances created by turbulent motion inside a jet 
ultimately radiates as sound. Since both radiating and 
non-radiating disturbances satisfy the same conservation 
laws, Lighthill’s equation is unable to separate the two. 
Additional criteria are necessary, such as a frequency- 
wave number, co-K, analysis of the turbulent motion and 
imposition of a condition of supersonic convective speed 
with respect to the ambient sound speed: co/k >a 0 . The 
traditional auto-correlation function in equation 4 is 
created by correlating two source points. In order to 
determine which part of this correlation ultimately 
radiates as sound, an additional wave number - frequency 
analysis of the right hand side is required (Morris et al 
2002, Goldstein, 1976). On the other hand, multiplication 
by the field point in equation 5, effectively imposes a 
filter function, since by definition a microphone kept in a 
far-field location only senses the radiated part of 
disturbances. Therefore, the cross-correlation function 


Rpuu; P - separates out hydrodynamic and acoustic 
fluctuations. There is no need to perform additional 
frequency-wave number analysis. 

The causality relationship requires a second time- 
derivative, which ought to be avoided in experiments. 
Fourier transform of equation 5 yields (noting that the 
autocorrelation function is transformed to power spectral 
density and cross-correlation to cross-spectrum): 


7tf df 


v puu.p' 


'(X f ,X s ;f) d 3 X s 


cross- spectrum : S puu p - (X f ,X S ;f) = (6) 


L R PUU ; P ' (X, , X f ; x) exp(- j 2 7tf x) dx 
Finally, the acoustic intensity at the field point X f : 

I(X f ) = - [ S , 2 df 

« 0 Po -i p 


71 

r Po a o 


Jjf 2 S„ p -(X„X,;f)d ! X,df 

-ooV 


(7) 
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The above equation shows that the intensity of sound 
radiation is directly dependent on the distribution of 
cross-spectrum function in the flow field. 

There are multiple issues surrounding the 
causality method. At a first glance the right hand sides of 
equations 4 & 5 are expected to be equal although the 
integrands are different. However, they may not be so as 
the difference between propagating and non-propagating 
disturbances is included in the latter. There is an issue of 
non-uniqueness in the application of causality method 
(Ffowcs Williams 1973). Since the sound pressure at the 
field point is a large sum over the entire sound-producing 
region of the jet, an unlimited number of variations in the 
source correlation can lead to the same summation at the 
field point. This criticism is not special for the causality 
technique but in general true for the more common form 
of source description via two-source points correlation 
(eqn. 4), which likewise, has to be integrated over the 
source volume. In a broader sense, many inverse 
problems in physics are of this nature. Although 
mathematically this appears to be a problem, it can be 
argued that the distribution of correlation functions 
measured in a real experiment is the valid distribution. 

The interpretation of the correlation function is 
another issue that brings back a long-standing discussion 
of whether turbulence consists of uncorrelated small 
eddies or large vortical structures of significant spatial 
coherence. If turbulence is described by the former then 
the <puu,p > correlation function measured from various 
points in the jet are mutually independent; the integration 
over the jet volume becomes a simple addition in power 
and the correlation data can be used to determine source 
efficiency along the jet axis (Seiner & Reethof, 1974). 
Such an interpretation also leads to a discussion of the 
number of eddies responsible for sound generation at a 
given instant (Lee & Ribner, 1972). It is now well 
established from various experimental observations that 
low Strouhal frequency turbulent fluctuations are 
dominated by organized vertical waves with significantly 
long spatial coherence. The effect of source-coherence on 
jet noise has been discussed by Michalke (1983). For 
sources with long coherence length scales the phase of 
<puu,p > correlation will vary from point to point and the 
integration over the jet volume requires a knowledge of 
both magnitude and phase. 

So far discussions have been confined to the 
contribution from the puu term only, ft is straightforward 
to extend this analysis to include the second term of Ty 
tensor, equation 1. The far field pressure fluctuations 
caused by pressure and density fluctuations in the plume 
are: 


p'(X f ;t) 

-Vl^Rp'-""VK- t '>] d3x . 

: a n r i d t 


( 8 ) 


471 i 


To completely evaluate this term one needs to measure at 
least temperature and density fluctuations inside the 
plume. The present experimental setup is unable to do the 
former. Neglecting contribution from entropy 
fluctuations, local sound speed a s relates pressure 
fluctuations to density fluctuations: p’ = a s 2 p’. Since a s 2 / 
do 2 ~ Tj/To, where T, is jet temperature and T 0 is ambient 
temperature, the above equation becomes: 
p(X f ;t) 


1 r d 2 

47t r * 3 1 /2 


U--1 P (X.;0 


d 3 X„ 


(9) 


where T 0 is the ambient temperature. For an unheated low 
Mach number jet T s /T 0 ~ 1 and the right hand side has 
negligible contribution to the far field. However, as the 
Mach number M increases the effectiveness of this term 
is expected to increase as M 2 . A reasonable estimate may 
be obtained by using time average plume temperature: 


M 


-1 


p (X f ;t) = L f A^[p' (x ;t ')] d 3 X s . (10) 

47t r J dt 

V 

Once again, multiplying the source integral with the far 
field pressure fluctuations: 

p' 2 (X,;t)) 

f . 


TJ 

vT 0 


-1 


4 *wM R A x .< x ^ )ld ’ x . <“) 

where, R p;p -(X s ,X f ;x) 

= (p / (X s ;t + x-T 0 )p / (X f ;t)) 

Since data analysis is performed in narrowband 
frequencies, Fourier transform of above yields: 


S p ,2 (X f ;f) df 


A 


1-^ 


7tf df 


0 J 


J v S PiP '(X f ,X s ;f) d 3 X s 


( 12 ) 


cross- spectrum : S <(X f ,X s ;f) 


= L R p-p' ( Xs ’ X f ; T ) ex P( - J 2 n f ' 0 dx 

Therefore, a correlation between the flow density 
fluctuations and the far field sound pressure fluctuations 
provides an estimate of contribution from the second term 
of LighthilFs stress tensor. 
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It is best to present the experimental data in a 
non-dimensionalized form. Following the traditional 
acoustic analogy approach, two different velocity scales 
are applied: ambient sound speed ao for field points and 
Uj for source points. The length and time scales are 
derived using the jet diameter D. Ambient value p 0 is 
used to normalize density. Using superscript * for non- 
dimensionalized quantities, S t for Strouhal number, and 
M ;i for Mach number based on ambient sound speed, 
equation 6 becomes: 


S*, 


,(x; ; s t )ds t 


-M^ds,J v s; unp .(x;,x; ; s,) d 3 x; , 


s.=®. 

U 


M„ - 


U; 


s ' T 


s ,- 


(13) 


puu, p 


puu,p , TT 2\/ 2 \ 

(PoUj ) (po^O ) 

Similarly, non-dimensionalized contribution from the 
second term (equation 12) is: 

s* (x; ; s t )ds t 


i-i 

T, 


M ; 


71 S 


S r,r' = 


o J 

s 


-^ds t { v s; p ,(x;,x; ; s t )d 3 x;,(i4) 


p,p 


(Po«o) 

In the present experimental program the principle idea of 
correlating flow fluctuations to the sound pressure 
fluctuations has been extended to include some other 
variables, which do not necessarily abide by the 
framework of Lighthill’s equation. For example, efforts 
are made to correlate pvv fluctuations not only with 
microphone signal from 90° to jet axis but to 30° position 
as well. Inquires are made as how do various other 
quantities, such as various Reynolds decomposed terms 
of puu: puu^ , p'uu" etc. correlate with sound 
pressure fluctuations. These additional efforts resulted in 
some interesting results. 


Fundamentals of flow measurement using Rayleigh 
scattering principle: 

A simplistic description of the measurement process, 
using laser induced Rayleigh scattering is schematically 
shown in figure 2. When a laser beam is allowed to pass 
through a gas, the molecules present in the gas cause 
inelastic and elastic light scattering. The inelastic part is 
called Raman scattering and the elastic part Rayleigh 
scattering. The Rayleigh scattering process describes most 
(~ 99%) of the molecular scattered light. In the present 
experiment scattered light is collected and spectrally 


resolved to measure velocity. Since the Doppler shift 
frequency is relatively small, a narrow linewidth incident 
laser beam is necessary to resolve Rayleigh spectrum. 
Even if the gas medium is stationary, the random thermal 
motion of the gas molecules creates a wide range of 
Doppler shift - resulting in a spectral broadening of the 
collected light. The Full Width at Half Maximum 
(FWHM) depends on the distribution of molecular 
velocities and, therefore, is a measure of gas temperature. 
In the case of a moving gas media, the bulk motion is 
superimposed on the random velocity of the individual 
molecules; therefore separation between the peaks of the 
incident laser line and the Rayleigh spectrum provides a 
measure of the bulk velocity. Finally, the total light under 
the Rayleigh spectrum is proportional to the molecular 
number density and provides a measure of gas density. A 
single Rayleigh spectrum carries information of one 
component of gas bulk velocity, temperature and density. 

This basic principle has been used in the past to 
measure time-averaged quantities (Panda & Seasholtz 
1999, Forkey, Lempert & Miles 1998, Elliott & 

Sammimy 1996). Since density variation modulates the 
total scattered light, unsteady density fluctuations are 
easier to measure (Panda & Seasholtz, 2002 among 
others). Extension of the Rayleigh scattering technique to 
measure unsteady velocity fluctuation has remained a 
challenge. Recently, Seasholtz, Panda & Elam (2002) 
have simultaneously measured the time variation of 
density and velocity fluctuations in free jets. Earlier 
efforts that culminated towards the present setup can be 
found in Seasholtz, Panda & Elam (2001) and Seasholtz 
& Panda (1999a, 2000). The present setup is for a point- 
measuring system; a continuous wave laser was used, and 
scattered light from a point on the beam was collected 
and analyzed. (It was realized early on that the current 
limitations of tunable pulsed- lasers and camera systems 
make an area measurement technique unsuitable for 
dynamic measurements over a large frequency 
bandwidth.) The scattered light collected from the probe 
volume was spectrally resolved by a Fabry-Perot 
interferometer. To illustrate the approach, first the nature 
of a fringe formed by the interferometer is shown in 
figure 3. The field of view in the image covers a fraction 
of the free spectral range. At first, a small portion of light 
directly out of the laser beam is imaged through the 
interferometer. The narrow line width of the laser makes 
a sharp, narrow ring in this image (reference image, 
figure 3a). In contrast, when the Rayleigh scattered light 
from a moving gas medium is imaged (Rayleigh image, 
figure 3b), a different, and diffused ring results. The 
radial shift in the peak intensity location between the 
Rayleigh and reference images corresponds to the 
Doppler shift associated with the bulk motion of the air 
stream (the diffused nature of the Rayleigh image is due 
to the thermal broadening). Since the laser frequency, and 


NAS A/TM— 2004-2 12385 


5 



therefore, the reference fringe is fixed, any radial shift of 
the Rayleigh image indicates an instantaneous change of 
air velocity. To continuously monitor the ring diameter in 
the Rayleigh image, Seasholtz et. al. (2001, 2002) used 
an image dissector that split the image into two 
concentric parts (figure 4), and measured the ratio of light 
intensity from the inner and the outer parts using two 
photomultiplier tubes (PMT). A calibration of the 
intensity ratio from known jet velocities was necessary 
for later use in unknown flows. To measure density and 
velocity simultaneously, a total of 3 PMT were used: 2 
for measuring velocity as described above and a third one 
to monitor scattered light power variation which is 
proportional to the air density fluctuation. 

The component of velocity measured using a 
given optical arrangement depends on the angular 
position of the collection lenses with respect to the 
incident laser beam. Figure 5 shows 2 different optical 
arrangements of the transmission and collection optics for 
measuring u and v velocity components. The Rayleigh 
scattered light was collected and transmitted to an 
adjoining room for spectral analysis. Figure 6 shows a 
schematic of the spectroscopic arrangement which was 
identical for both velocity components. Since density and 
a component of velocity was measured simultaneously, 
either p(t), u(t) or p(t), v(t) were measured using the two 
different collection arrangements. Figure 5(b) presents a 
scattering diagram pertinent to the radial v component 
measurement, k, is the incident wave vector normal to the 
jet flow direction, k is the scattered wave vector pointing 
towards the collection optics and, k = ks - k ; , is the 
scattering vector. The arrangement measures Doppler 
shift, f d = kU/27t from the radial velocity component v. 
Similar diagram for the optical arrangement of figure 5(c) 
shows that the measured Doppler shift corresponds to -u 
component. The calibration process accounted for the 
sign reversal. 

The image formed by the Fabry-Perot 
interferometer depends on the optical frequency 
distribution of the scattered light, properties of the 
interferometer, called instrument function, diameter of 
the optical fiber used for transmitting collected light and 
focal lengths of the collimating and fringe forming lenses 
of figure 6. A detailed discussion can be found in 
Seasholtz, Panda & Elam (2002). In brief, the optical 
spectrum of the Rayleigh scattered light 91(v-Vo) carries 
information on the velocity distribution of gas molecules. 
For a low density gas the spectrum is given as: 

9t(v- v 0 )df 



2n(v- v 0 )-k.U 

a P k 



( 15 ) 


Where, Vo is the incident laser frequency and a p , the 
“most probable molecular speed,” is related to the local 
sound speed a and ratio of specific heat y through 
a p =(2/y) 12 a. An examination of the Gaussian distribution 
shows that the width of the spectrum depends on a p k, 
which is related to Vt and the peak is shifted by the 
Doppler shift kU. The image formed, after the Rayleigh 
light is passed through a Fabry-Perot interferometer, is 
basically a convolution of the Rayleigh spectrum with the 
instrument function l F p(v, 0 r ). The light power 
distribution Pi at any position (r, 0) in the image plane is 
given as: 

P : (r, 0) = P^k J J SR(v) !fp ( v A- ) dv dA (16) 

Where, Pj; is the total Rayleigh scattered power, and 0 r is 
the angle made by a light ray reaching the elementary 
area dA in the image plane with optical axis. In the 
current setup this image, formed at the focal plane of the 
fringe-forming lens, was dissected into two parts by a 
concentric, tilted mirror assembly (image dissector in 
figure 6) and measured by two PMT. The intensity 
variations from the two parts were measured from known 
velocity flows. It was found that the intensity variations 
in either part of the image can be modeled by 2 nd order 
polynomials. If N 2 and N 3 denote photo-electron count 
rates from the inner and outer PMT then: 

N 2 = zfj +u B t +u C ; , N 3 = A o +u B o +u C 0 (17) 

Where, A;, B ; , Q, A,, B 0 and C G are calibration constants. 

The velocity component u is measured from a ratio of the 
two counts, R=N 3 /N 2 : 

_ - (RBj -B 0 ) + 4{RB, -B o f -4 {RA, -A o ){RC, -Cj 
2(RCj-C 0 ) 

( 18 ) 

Note that velocity is calculated from a ratio of two counts 
R=N 3 /N 2 , which cancels out changes in overall scattering 
intensity associated with flow density variation. Also the 
physically meaningful, negative root of the quadratic 
equation is considered. The above analysis does not 
account for changes due to temperature fluctuations on 
Rayleigh spectra. A numerical uncertainty analysis 
(Seasholtz et. al. 2001), however, demonstrates that such 
an effect is small compared to the bigger change 
associated with velocity fluctuations. 

Density variations were measured by splitting off 
about 1 0% of the collected light using a beam splitter 
(BSl in figure 6) before the rest is passed through the 
interferometer. The power variation in this split part is 
measured by PMT1. It is known that for a fixed optical 
setup and a fixed composition of gases the total scattered 
light is directly proportional to the gas density: P 3 ; °< p. If 
the photo-electron count rate from PMT 1 is A then: 

N 1 = Cpl + PCp2 > ( 19 ) 
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Where, C p[ and C P 2 are calibration constants determined 
from measurements in know density flows. 

II. EXPERIMENTAL SET-UP 

Experiments were performed at NASA Glenn Research 
Center using three different nozzles (one convergent and 
two convergent-divergent) operated at Mach numbers, M 
= 0.95, 1.4 and 1.8. The convergent-divergent nozzles 
were designed by the method of characteristics and their 
geometries were reported in Panda & Seasholtz (1999). 
The operating conditions are described in Table 1. All 
nozzles were 25.4 mm in exit diameter. The jet facility 
used a continuous supply of unheated compressed air. 

The facility was located in a large test chamber, which 
was not anechoic per se, but acoustic absorbent material 
was placed around the vicinity of the nozzle and in the 
ceiling and walls of the test cell to minimize reflection. 
Two 1 /4-inch microphones were used to measure sound 
pressure fluctuation spectra. The microphones were 
mounted on an arc that allowed positioning at a distance 
of 50 diameters from nozzle lip and polar angles from 30° 
to 90° to the jet axis with 10° increments. For the bulk of 
the experiment one of the microphones was kept fixed at 
30° to the jet flow direction and the other at 90°. The 
Rayleigh scattering system is somewhat elaborate and the 
following provides a description. 

An in-depth discussion of the Rayleigh set-up 
can be found in Seasholtz, Panda & Elam (2001, 2002). 
Following is a brief discussion of some of the central 
features. The optical system was built in two parts. The 
first one is around the jet facility for transmitting laser 
light and collecting the scattered light (figure 5). The 
scattered light was then passed through 0.55mm diameter 
optical fiber to a quiet room where the second part, 
consisting of spectroscopic and photon counting 
electronics, were placed (figure 6). The splitting of the 
setup is a part of special attention that had to be taken for 
successful implementation of the optical setup. Care had 
to be taken to reduce dust particles in the air streams, to 
minimize the effect of vibration on the optical equipment 
and to stabilize the interferometer from temperature and 
vibration induced drifts. To reduce dust particles the dry 
air, supplied to the facility from a central high-pressure 
facility, was passed through sub-micron filters. Thi s 
made the primary jet air very clean. To clean the 
entrained ambient air, an additional air blower & filter 
system was installed that provided a large, 200mm 
diameter, low-speed (~ 20m/s) co-flow around 25.4mm 
primary jets. 

There were two different optical systems built 
around the jet facility to measure u and v components. 
Both of these were built over an X-Y traversing unit that 
carried laser head, transmission and collection optics. The 
probe volume was moved from point to point in the 
plume, and thereby, allowing survey over a cross- 


sectional plane. To maximize utilization of the available 
light the laser head of the solid-state, frequency-doubled, 
Nd:VC >4 laser was placed at the bottom part of the set-up 
(figure 5a). About 5 watts of single mode, 532nm 
wavelength laser light was transmitted through a hollow 
side beam that contained a half-wave plate, focusing lens, 
mirrors and baffles. Since Rayleigh scattered light is 
polarization dependent, the half-wave plate was rotated to 
align the peak scattering plane with receiving optics. The 
background scattered light was significantly attenuated 
by suitable use of baffles and beam-dump. The laser head 
was placed in the same traversing unit as the rest of the 
transmission and collection optics. It was found that the 
noise from the jet created a tonal excitement of the laser 
line at around 430Hz. To reduce this excitation an 
anechoic box was built around the laser head. This box 
significantly reduced the laser unsteadiness, but a trace 
remained and manifested in the experimental results. The 
incident laser paths and positions of the collection optics 
differed for the two setups to measure radial (figure 5a) 
and axial (figure 5c, 5d) velocity component. For the 
former, laser beam was passed normal to the jet axis and 
the collection optics were placed vertically down to 
collect scattered light from 90° to the incident beam as 
well as jet axis. On the other hand, to measure axial 
velocity component the laser beam was passed at 45° to 
the jet axis, and 90° scattered light was collected by 
lenses oriented in the same plane containing the laser 
beam and the jet axis. The lens combinations for the 
collection optics were the same for both setups. The 
scattered light was first collimated by a 300mm focal 
length, 80mm diameter achromat, and then focused by a 
160nrm focal length achromat on the face of a ,55mm 
diameter multimode fiber. The combination of the fiber 
diameter and the magnification ratio of the collection 
optics fixed the probe volume length to 1 .03mm. Not 
shown in figure 5 is an additional part of the set-up where 
a small part of light from the transmitting beam was split 
off for the purpose of monitoring laser frequency as well 
as to maintain a close-alignment in Fabry-Perot 
interferometer. Whenever necessary, a pneumatically 
actuated mirror was placed in the transmission laser path 
to divert light towards a diffuser. The diffuser was 
simultaneously brought in front of the collection fiber and 
reference light, instead of the Rayleigh scattered light, 
was passed through the collection fiber. 

To protect the spectroscopic equipment and the 
split mirror assembly from high noise level generated by 
supersonic jets, they were separately kept in an adjoining 
quiet room. This second part of the setup is schematically 
shown in figure 6. Fight arriving via optical fiber was 
collimated by 1 OOmrn focal length lens F 1 , and about 
10% was split by BS1 and measured by PMT1. Output 
from PMT 1 provided information on the modulation of 
total scattered light intensity and therefore, was 
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proportional to the air density fluctuations. The rest of the 
collimated beam was passed through a 70mm aperture 
Fabry-Perot interferometer for spectral analysis. Single 
wavelength light from an extended source is imaged as 
consecutive rings (fringes) at the output of the 
interferometer. However, restriction of field of view, 
imposed by the fiber diameter, created only one fringe as 
shown earlier in figure 3. The fringe-forming lens, which 
ultimately images the fiber face on the image dissector, 
was made of two camera lenses with suitable separation 
for an effective focal length of 2909mm. The large 
magnification ratio of the setup created a 16mm diameter 
image of the fiber on image dissector. The image 
dissector was made of two concentric, and slightly tilted 
mirrors. The inner one had a diameter of 10mm and 
directed the inner part of the fringe to PMT2, while the 
25mm diameter outer one directed the outer part of the 
fringe to PMT3. The ratio of light intensities from PMT2 
and PMT3 provided a measure of either axial or radial 
velocity as described earlier. 

The success of the velocity measurement system 
was critically dependent on stable operation of the 
interferometer. Slight thermal drift or change in the 
incident laser frequency displaced the reference fringe 
(figure 3a) and manifested as an artificial bias in velocity 
measurement. This made the auto-alignment setup (figure 
6) a necessary part. The alignment system was basically a 
feedback control that first measured the reference fringe 
diameter and compared it with a prescribed targeted 
diameter. This function was accomplished by splitting 
parts of the transmitted light using a 3-prism assembly, 
and imaging them on a CCD camera. Subsequently, the 
difference between the targeted to the measured fringe 
diameter was translated into a change of high voltage 
supply to the piezo-electric actuators that adjusted 
interferometer’s plate separation. Before every Rayleigh 
measurement, reference light was collected, passed 
through the same fiber and optical system, and the auto- 
alignment system was engaged. When the desired fringe 
diameter was obtained within a tolerance, the reference 
light collection system around the jet facility and the 
prism assembly were disengaged, and velocity and 
density measurement via analysis of Rayleigh scattered 
light began. 

Photoelectron counting electronics were used 
with all 3 PMT signals. The counting was performed over 
contiguous bins of prescribed time-duration. The 
digitization of the microphone signals was synchronized 
and finally 3 channels of photoelectron count and 1 
channel of microphone signals were simultaneously 
acquired in a Personal Computer. Photo-electron 
counting was started at the beginning of a clock pulse that 
also digitized the analog microphone signal. The counting 
ended at the next clock pulse. Long time records of up to 
5 million data points from each PMT and the microphone 


channel were acquired with a typical sampling rate of 
90,000/sec. The data points were converted to physical 
parameters: instantaneous density, velocity and sound 
pressure through the usage of proper calibration 
constants. The Welch method of modified Periodograms 
(1967) was used to calculate individual spectrum and cross- 
spectral density. Each long record was divided into small, 
50% overlapped, segments; modified periodograms of each 
segment provided local estimates; average of all local 
estimates provided the final cross-spectral density. The 
correlation calculations were also performed via Fourier 
transform, where individual segments of velocity and sound 
pressure time histories were Fourier transformed, 
multiplied and inverse transformed. Average of local 
estimates produced the final correlation. 


III. RESULTS: 

Calibration and validation: As described before, 
calibration constants are required to convert photo- 
electron counts from the 3 PMTs to density and velocity 
values. For this purpose measurements had to be made in 
known velocity and density flows. The axial velocity and 
density in the potential core of jet plumes from a 
convergent nozzle operated in the subsonic Mach 
numbers, and CD nozzles at the correct operating 
conditions, are known from isentropic relations: 


RTo 


Pay 


; u = Ma 


where, M = 


Y 


.PaJ 


y— 1 A 


and a = J yRT 0 
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( 20 ) 
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The ambient pressure P a , plenum pressure Po and Plenum 
temperature To were monitored using pressure- 
transducers and a thermocouple; for unheated air y=1.4. 
Since the jets exhausted into ambient, the plenum 
pressure was changed to vary plume Mach number. Axial 
velocity calibrations were performed in the same plumes 
where later measurements were performed. For the 
calibration of the radial component of velocity, however, 
a separate small nozzle was mounted close to the probe 
volume. This calibration jet was rotated to align the 
plume along the radial velocity direction of the primary 
jet. Figure 7 shows a set of typical calibration curves for 
density and axial velocity measurements. Parts (a) & (b) 
show variation of count rate from the two PMT that 
measure light intensities from the split images formed 
after the Interferometer (PMT2 and PMT3 in figure6). 
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The solid lines are least square fit of a 2 nd order 
polynomial as described in equation 17. Figure 7(c) 
shows ratios of the two counts, measured at different flow 
velocities along with a solid line representing calculated 
velocities (equation 18). This plot shows a reasonably 
good fit to the experimental data. Finally, density 
calibration is shown in figure 7(d). One noticeable feature 
is a relatively high scatter of the calibration data from the 
supersonic Mach 1.4 and 1.8 plumes. Prior density 
measurements and schlieren photographs in these plumes 
(Panda & Seasholtz, 2002) showed the presence of weak 
shock structures even at the design operating condition. 
Perhaps weak shocks are impossible to avoid in any free 
supersonic shear flows. To average out the density and 
velocity variations across the periodic shock system, data 
was obtained from various axial stations within the 
potential core. The data scatter of figure 7 is a reflection 
of point to point variation in the supersonic plume. 

Naturally, this data scatter has added some 
uncertainty in determining turbulent properties. There 
were additional minor sources; however, the combined 
effect of all such sources was masked by the uncertainty 
from electronic shot noise, inevitable in any optical 
intensity measurement. This is illustrated in figure 8 
where velocity fluctuations spectra measured using the 
present Rayleigh arrangement is compared with that 
obtained using a hot-wire probe. This figure shows that 
the shapes of velocity spectra are similar while the 
absolute energy level in the spectrum measured by the 
Rayleigh technique is two and a half times of that 
measured by hot-wire. The Rayleigh spectra float on 
constant white noise floors, which are expected 
consequences of electronic shot noise. The propagation of 
electronic shot noise in the velocity spectrum is 
somewhat complicated by the ratio of photoelectron 
counts, R=N 3 /N 2 , needed to determine instantaneous 
velocity. The shot noise randomly changes the ratio and 
manifests as a noise floor in the velocity spectrum. This 
constant floor is all that exists in data taken from still 
ambient air. It should be mentioned that the primary 
interest of the present work is in cross-correlation 
between sound pressure fluctuations and turbulent flow 
fluctuations. Since the noise sources in the two 
measurements are different, the cross-correlation process 
eliminates the influence of electronic shot noise and the 
absolute cross-correlation values are expected to have 
much smaller level of uncertainty. 

Nevertheless, a reasonable estimate of the 
standard deviation (root-mean-square) of various 
fluctuating quantities is required to non-dimensionalize 
the absolute correlation values. For example, velocity 
fluctuations to far field sound pressure fluctuations are 


expressed as: 


( U ’P / ) 

^rmsP rms 


, which requires an estimate of 


Urms. Such estimates were calculated from spectral data by 
subtracting the constant shot noise floor shown in figure 
7. The noise floor was estimated as the spectral density 
value at the highest resolved frequency of 50KHz, 
(P u )f=5oooo, where it is known that the energy from 
turbulence fluctuations is small. 

= V Af Z( P . PI) 

Flere, P u represents power spectral density of u 
fluctuations measured in (m/s) 2 /Hz and Af = 50,000 is the 
frequency range. The rms values measured by this shot 
noise subtraction process are found to be reasonable. For 
example, in figure 7(a) (u nns /Uj) ho t-wire = 0.166 and 
(u rms /Uj) R a y icigh corrected = 0. 154 and in figure 7(b) 

Ofrms/Fij)hot -wire 0.14 and (u 

rms/Uj)Rayleigh corrected 0.132. 

This subtraction method was used uniformly to estimate 
Prms, Urms and (puu) m]s from respective spectra. 

Density and velocity fluctuations spectra: Figure 9 show 
typical velocity and density spectra obtained from the 
strongest noise-generating regions of Mach 1.4 plume. 

All velocity data is non-dimensionalized by the centerline 
axial velocity Uj while density is non-dimensionalized by 
the ambient value. The puu spectrum of figure 9(d) was 
obtained by first multiplying instantaneous density and 
velocity: puu(t)=p(t)u(t)u(t) and then taking Fourier 
transforms of the resulting time series data (note that 
mean values p and u were not subtracted). The density 
velocity cross-spectrum was obtained by multiplying 
Fourier transform of individual time traces: 

S up (f) = 3 (p(t) ) * 3 * (u(t) )= Re{S up }-i Im{S up } 

magnitude |S up | = ^(Rcfs^}) 2 +(lm{S up }) 2 , (22) 

Im{s } 

phase9 "' = ^l 

In this equation complex conjugate is represented by the 
superscript *. The magnitude of the cross-spectral density 
is plotted in figure 9(e), and phase in 9(f). Note that the 
sharp spike at very low Strouhal frequency (430Hz in 
absolute frequency) is due to a spurious oscillation in the 
laser frequency and should be ignored for further 
interpretation. It is possible to make multiple interesting 
observations from figure 9. First, density fluctuation 
spectrum is similar to that of axial velocity, while radial 
velocity shows some difference: u and p spectra show a 
continuously decaying shape, while v-spectrum has a 
hump shape. This indicates a difference in eddies that 
produce most energetic fluctuations in axial and radial 
velocity components. However, data obtained from the 
shear layer show similarity in shapes for all 3 spectra: p, 
u & v. That the density fluctuation spectra follows that of 
axial velocity is somewhat expected, as u fluctuations are 
more energetic than v fluctuations. The difference 
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between u and v spectra, however, is unexpected. In the 
normal mode representation of turbulent fluctuations, 
which is commonly employed in hydrodynamic stability 
calculations, it is customary to assume identical spectral 
distribution for all flow and thermodynamic variables. 
This obviously is incorrect, at least around the centerline 
region of the jet, which the present experiment finds to be 
the strongest noise source. A second observation from 
figure 9 is that the puu spectrum, as expected, has the 
same shape as that of u spectrum although in absolute 
value (puujrms » u nns . Third, the velocity-density cross- 
spectrum has a continuously decaying shape, indicating u 
and p fluctuations are the best correlated at the larger, 
lower frequency eddies and progressively de-correlate as 
the eddy size becomes smaller. Fourth, the flat phase 
relation of figure 9(f) implies that density and axial 
velocity fluctuations are in phase for all eddy sizes. 
Following is a discussion of measurement uncertainties. 

As already mentioned, the fundamental source of 
uncertainty in optical measurement is due to electronic 
shot noise. Since density spectra are calculated directly 
from the Fourier transform of photoelectron counts, the 
effect of electronic shot noise is determinable. A detailed 
discussion can be found in Panda & Seasholtz (2002). In 
essence, like velocity spectra, shot noise adds a constant 
floor and randomness in density spectra. The latter is 
significantly reduced in the present experiment by 
averaging over many segments obtained from 5 million 
point long data string. The constant noise floor 
particularly masks the lower energy, high-frequency side 
of the spectrum. Panda and Seasholtz (2002) used a 2- 
PMT-correlation technique to alleviate this problem. In 
the present work no such correlation technique is applied. 
It is estimated that as much as 50% of spectral energy at 
St = 0.2 may be due to this electronic noise source. 

Aside from the shot noise, there are three other 
smaller sources of noise contributing to velocity spectra. 
First, a slow random variation in the laser frequency over 
30 MHz (.001 cm' 1 ) that translates into ±8m/s velocity 
fluctuation. Second, the intense sound level produced by 
the jet caused a periodic variation of the laser frequency 
at about 430Hz. As described earlier the laser head was 
enclosed in an anechoic box to alleviate this problem. 
Nevertheless, the periodic variation could not be 
completely eliminated; as evident from u, v, puu spectra. 
Finally, the Fabry-Perot stabilization setup, needed to 
lock the reference fringe at a fixed radial position, was 
effective only within an error margin. A change in the 
reference fringe position, from the value used for 
calibration, translates into a spurious mean velocity in the 
analyzed signal. Similarly, a random positioning error in 
the reference fringe results in added energy in the spectral 
data. A direct estimate of the uncertainty from all such 
sources is difficult to obtain. Nevertheless, the success of 
the causality method is critically hinged on noise 


cancellation obtained in cross-correlating two signals of 
independent noise sources. This cross-correlation is 
described in the next section. 

Correlation between flow fluctuations and sound 
pressure fluctuations: The presentation of correlation 
data starts with figure 10, where normalized cross- 
correlation plots for 30° and 90° microphone angles are 
presented. The cross-correlation values were calculated 
via Fourier transform, that is, the cross-spectrum was 
calculated first, and then an inverse transform was taken 
to return to the time domain. Neither the flow data, nor 
the microphone data was shifted in time. The cross- 
correlation data show sharp rise at a time delay x 0 , which 
is found to be the time needed for sound waves to travel 
from the laser probe location to the microphone location, 
To= r/ao. This confirms that the “cause” (p, v, pvv 
fluctuations in the jet) produces the “effect” (microphone 
pressure fluctuations) at a predictable time lag, and in 
turn provides confidence on the measured data. Figure 10 
presents data from 2 different Rayleigh setups: parts (a), 
(b) & (c) are from u-setup, parts (d) & (e) from v-setup. 

In either case the laser probe was located at the end of the 
potential core, which was found to be the strongest sound 
source. There are multiple interesting observations that 
can be made from figure 10. First, sound pressure 
fluctuations p' correlates far better with u or puu 
fluctuations than v or pvv fluctuations. Second, air 
density fluctuations inside the jet show as good a 
correlation as puu. Third, sound pressure fluctuations at 
90° to the jet axis correlate very poorly with any flow 
variables. The <pvv;p > correlation of figure 10(e) is of 
particular interest. The second derivative of pvv 
represents a longitudinal quadrupole with peak radiation 
direction 90° to the jet axis. The radiation from such 
quadrupoles is unaffected by convective amplification. 
Therefore, it had been expected that <pw;p > correlation 
from 90° microphone would be significant. This 
obviously was found to be incorrect. Since, v and pvv 
fluctuations are found to be poorly correlated with far- 
field noise, such data are not presented in the rest of the 
paper. Fourth, the time duration Ax over which 
correlation changes from zero to negative to positive and 
back to zero is significantly long (Figure 10c). The 
trailing edge of the correlation is somewhat difficult to 
determine due to the spurious ripples arising from the 
reflected waves (discussed in the next paragraph); 
nevertheless for the particular jet, Ax is estimated as 
0.53milli-sec. This provides a measure of coherence time 
of the turbulent eddies responsible for sound radiation. 
Invoking Taylor’s hypothesis, the corresponding 
coherence length scale is l CO her = Ax U c , where U c is the 
convective speed. Assuming U c =0.7Uj the coherence 
length lcoher ~ 6D. In other words, the longest eddy 
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contributing towards correlation is, on an average, 6 jet 
diameters long. Therefore, the noise source responsible 
for sound radiation, at the shallow 30° angle, is indeed 
“non-compact.” Following is a discussion of uncertainty 
in correlation data. 

As mentioned earlier, the cross-correlation process 
significantly reduces the shot noise contribution and 
therefore absolute values of <p;p> and <puu;p > are 
relatively error free. The primary source of uncertainty is 
due to the shot noise elimination process used in 
estimating root-mean-square values: p ms , (puu) rm „ etc. 
The shot-noise subtraction process has been described 
earlier. From repeated measurements and comparison 
with hot-wire data it is estimated that there is a ±10% 
uncertainty in rms calculation, which produces an equal 
uncertainty in correlation data. A second source of 
uncertainty is sound reflection from large lenses placed 
close to the plume. Although most of the optical 
components and mounts were covered by polyurethane 
foam to minimize reflection, some could not be covered 
for obvious reasons. The reflected sound waves arrive at 
the microphone position at a longer time delay than those 
reaching directly from the source. The additional ripples 
after the large primary spike, seen in figure 10, are due to 
the delayed reflected waves. For the most cases, such as 
in figure 10, it is easy to separate out the spurious 
reflected part from the desired direct correlation and is 
expected to have minimum effect on the maximum 
correlation value. 

It is difficult to cross check the correlation values 
measured in the present experiment with many of the 
earlier experiments. As mentioned earlier, a large number 
of similar investigations in the 1970’s used intrusive 
probes, such as hot-wires or microphones. An 
examination of data from these experiments shows 
unreasonably large correlation (in very low Mach number 
jets), which is rightfully assigned to the additional noise 
created by probe intrusion. One of the exceptions is the 
work of Schaffer (1979), who used Laser Doppler 
Velocimetry to measure u-fluctuations in Mach 0.98 jet 
and cross-correlated with sound fluctuations measured by 
a microphone. The <u ;p > correlation reported by 
Schaffer is comparable to the similar, Mach 0.95, data 
obtained in the present experiment. For example, 

(u'lp 1 ) 1 

=0.03 , when the probe 

^ rms P rms 

- max 

was positioned at x=10D and centerline and microphone 
was located 30° to the axis. The same correlation number 
was obtained in the present experiment with identical 
probe and microphone locations. This close comparison 
has provided further confidence on the measured data. 


Reynolds decomposition of puu and correlation with 
sound pressure fluctuations: In high Mach number jets 
the largest term in the Lighthill stress tensor arises from 
puu fluctuations, and fig 1 0 confirms that indeed the 
highest correlation is measured from this quantity. From 
the modeling perspective it is of interest to break down 
the full stress term into time averaged and fluctuating 
parts and to see how the individual components correlate 
with the far field noise: 
p=p+p', u = u+u' 

(puir,p'^= (2puu'-p') + (p'uu; //)+ ( 23 ) 

(2up'u'; p'^ + ipii'ii'; p'^j + (p'u'u'; p'} 

The 5 fluctuating terms in puu decomposition were 
constructed from the measured p(t), u(t) time traces and 
individually correlated with the sound pressure 
fluctuations from the 30° microphone. Figure 1 1 presents 
this data along with correlation from the full puu term. 
Clearly, the largest contribution is from the “shear noise” 
term2 puu 7 , and the next important term is p'uu : the 
other 1 st order term. The terms which are 2 nd and 3 ld order 
in fluctuations contribute minimally. A closer 
examination of figure 1 1 reveals that the phase of the 
correlation is different between various order terms. The 

time delay at which the maximum in (2 p u u^p') 

correlation occurs is different from that in (2 p u , u , ;p / ) 

correlation. In general, correlation with first and third 
order terms are nearly in phase while the second order 
terms are opposite in phase. Therefore, a simple sum of 
the +ve maxima values in the constituent correlations do 
not add up to value measured in <puu;p > correlation. 

The relative magnitudes and phase relation in the 
correlation data is further explored in Table II, where 
numerical values of contribution from various terms at 
the time delay corresponding to <puu;p > ma x are presented 
for various Mach number conditions. The observations 
made in the earlier discussions are once again supported 
by the numerical values. 

From the jet noise modeling perspective, there 
are three significant observations, which may lead to 
further improvements of the existing models. The first 
observation is about the Lilley’s equation (1974) which 
has been used to model jet noise (Khavran, Bridges and 
Freund 2002, Goldstein 2001 and others). For a 
transversely sheared flow, i.e, a flow where all mean 
velocity and temperature gradients, except that along the 
radial direction du/dr, is absent, the Lilley’s equation can 
be written as (Goldstein, 1976, Colonius, Lele & Moin 
1997): 


Schaffer found 
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Dt dt dx y p 0 


In the above equation u,' = u', v', w ' are the fluctuating 
velocity components, x ; = x, r, x 3 is the co-ordinate 
system, p is static pressure which is normalized by a 
suitable reference p 0 . This is a convected wave equation 
with the right hand side representing the sound source. 
Note that the source terms are second order in 
fluctuations. Goldstein (2001) provides the significance 
of this arrangement: “The dominant part of the Lighthill 
source term is quadratic in the total flow velocity, which 
can be decomposed into a mean plus a fluctuating 
component. The source function therefore contains terms 
that are both linear and quadratic in the fluctuating 
velocity components. Lilley (1974) argued that the linear 
terms, which are typically much larger than the quadratic 
quantities, do not actually radiate any sound and should, 
therefore, not be included in the source function, since 
they would tend to dominate over the much smaller 
quadratic terms which are the true sources of sound.” The 
experimental data of figures 10, 1 1 and Table 11 clearly 
contradicts above conviction. Derivation of equation 24 
actually requires moving term linear in fluctuation to the 
left hand side of the equation and rearranging using the 
momentum equation. This step perhaps, has mixed the 
source terms with the propagation terms. There is a 
caveat to the above discussion: in the following it will be 
demonstrated that the correlation measured in the present 
experiment is primarily due to the large coherent 
structures; contributions from the fine grained turbulence 
fall below the experimental noise floor. Therefore, the 
above discussion of limitations strictly applies to the 
noise component produced from the large scales only. 

The second modeling issue is the frequently used 
assumption of constant flow density in evaluating the 
density*velocity*velocity product term, i.e., puu ~ puu. 
All earlier studies employing the causality technique also 
used the same assumption. The present data, in contrast, 
show that the neglected p'uu term contributes from 10% 
to 25% towards <puu;p'> correlation. Therefore, density 
fluctuations should be included in modeling of the 
turbulence stresses. The third and the final observation is 
related to the first and is about modeling the space time 
correlation <puu(X s ,T); puu(X s +^, t , +x)> required for 
calculating the far field sound (equation 4). Following the 
intuition mentioned earlier, models depend on the length 
and time scales calculated from the second order u'u' 


fluctuations: <u'u'(X s ,f); u'u'(X s +^, t'+x)> (for example, 
Khavaran, Krejsa & Kim, 1992). The present experiment, 
on the other hand, shows the largest contributions, 
towards the noise radiated close to the jet axis, are from 
the first order fluctuations u' and p'. Therefore, the 
required length and time scales should be calculated from 
<u(X s ,f ); u(X s +t t'+x)> and <p'(X s ,t'); p'(X s +^, t'+x)> 
correlations. This is true for predicting noise from large 
organized structures. 

Comparative study of correlation at various microphone 
angles: Figure 12 shows variations in the maximum 
positive correlation coefficient, when the microphone polar 
angle is varied from 30° to 90° while the laser probe 
volume was kept fixed at the indicated (x, r) location. The 
angles are measured from the downstream direction. The 
microphone radial distance from the nozzle was kept fixed 
at 51 diameters. Interestingly, sound pressure fluctuations 
from the shallowest polar angle of 30° show the highest 
correlation. The correlations drop sharply to the noise floor 
at 60° and higher angles. The present correlation study is 
identifying turbulent fluctuations that radiate close to the jet 
axis. The absence of correlation at higher polar angles also 
implies existence of a second mechanism of noise 
generation which the present point measurement technique 
is unable to detect. This is discussed at the end of this 
section. 

Comparative study of correlation from various parts of 
Jets: For this part of the study the microphone was kept 
fixed, at 50 diameters and 30° to the jet axis and the laser 
probe was moved from point to point in the plume. Figure 
13 shows the variation in peak correlation when the probe 
was moved along the jet centerline. The highest correlation 
values were measured just downstream of the end of the 
potential core, perhaps due to merging and fragmentation of 
large organized structures present in the peripheral shear 
layer. Note that the potential core lengthens with Mach 
number due to a decrease in shear layer growth rate. The 
end appears around x/D = 6,1 and 9 for Mach 0.95, 1.4 and 
1.8 jets respectively (Panda & Seasholtz, 2002). The 
variation of the correlation coefficient along the shear layer 
(r/D = 0.45) is very different from that along centerline. For 
figure 14, the microphone was again kept fixed at the same 
30° location, while the probe volume was moved axially 
from point to point along r/D = 0.45. These data show a 
shaip difference between subsonic and supersonic plumes. 
In the subsonic Mach 0.95 plume all correlations fall below 
experimental noise floor, while in the supersonic Mach 1 .4 
and 1.8 cases significant correlation is measured. The 
difference between the centerline and the shear layer 
behavior prompted a study of radial dependence of the 
correlation coefficients. Figure 15 shows that indeed the 
highest correlation is measured at the centerline. The 
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<puu;p> correlation progressively falls to the noise floor 
for r/D>0.6. The <p;p> correlation, however, shows an 
interesting trend: a low correlation, yet higher than the 
noise floor, is measured when the probe volume is 
positioned outside the turbulent plume. This small value 
remained nearly constant beyond r/D=0.6. A qualitative 
explanation is as follows. Sound waves can be thought as 
density fluctuations. In Mach 1.4 jet, where some turbulent 
eddies attain supersonic speed relative to the ambient, the 
shock waves associated with the individual eddies 
propagate directly to the far field. The presence of this 
Mach wave emission process, perhaps, results into the 
measurable correlation between the near field (existing just 
outside the jet boundary) air density fluctuations and the far 
field sound pressure fluctuations. 

Frequency analysis of cross-correlation data, cross- 
spectral density: So far all correlation data were presented 
in time domain. Figure 16 presents data in frequency 
domain. The cross-spectral density values are normalized 
by time averaged plume and ambient properties: do, Uj, and 
Po, as described in equations 12 and 13. Note that figures 
16(bl) and (b2) were obtained from the same data set used 
in the spectral plots of figure 9 and correlation plots of 
figures 10 and 11. The first point to be discussed is that 
figure 16 provides evidence to the claim that cross- 
correlation data are significantly free of the shot-noise 
effects. Unlike p and puu spectra in figure 9, the S p , p - and 
Sp U u,p' cross-spectra rise 2 to 3 decades above the noise 
floor. Moreover, the shaip peak at 430Hz seen in puu 
spectra is absent in S puu , p - cross-spectra. The second point to 
be made is that the turbulent fluctuations providing most of 
the correlations lie in the nominal Strouhal number range 
0< St <1.0, with the peak value around St = 0.2. All of 
these cross-spectra demonstrate that the measured 
correlations are mostly due to the large organized structures 
present in the plume. Contribution from small scale 
fluctuations fall below the noise floor. Another noticeable 
trend is the similarity in the overall shapes of S p , p - and S puu , p - 
cross-spectra over the entire Mach number range. 

Contribution to far field noise from <p;p> and <puu;p> 
correlation: The causality principle, described in the 
introduction, can be used to back calculate far field sound 
pressure fluctuations from the measured cross- 
correlations. This obviously is not a modeling approach, 
yet serves the purpose of estimating the fraction of noise 
generated from various regions of the jet. Figure 17 
presents one such effort to calculate auto-spectrum of 
sound pressure fluctuations at the 30° location (and 50D 
away from the nozzle exit) from a single point cross- 
correlation measurement. The laser probe was placed at 
the centerline and downstream of the potential core where 
maximum correlation is measured. The probe and 


microphone locations are the same for which correlation 
and spectral data are presented earlier. Figures 17(c) and 
(d) show a comparison of the calculated auto spectrum 
with the actual measured profile. Recall that the 
unevenness in the measured auto-spectra is an artifact of 
reflection from large lenses placed close to the plume. 

The calculation procedures involve multiplication of 
cross spectral densities, S p , p - and S puu , p , by various terms, 
most importantly by frequency squared, as described in 
equations 13 and 14. The dashed curve represents left 
hand sides of the equations and the solid curve the right 
hand side. Note that the causality equations require an 
integration of all cross-spectra measured over the entire 
jet plume. This has not been attempted in the present 
work. Instead the calculated auto spectrum should be 
interpreted as effectiveness of a unit volume of flow, 
centered at the laser probe, in creating sound pressure 
fluctuations at the microphone location. Interestingly, 
figure 17(d) shows that the measured correlations are so 
high that such a unit volume by itself produces more 
noise than the measured auto-spectrum. Figure 17(c) 
shows that the dipole contribution from density 
fluctuations is about 1 0% of the total measured sound 
pressure. Perhaps the only way the extraordinary 
efficiency of a single volume of turbulent fluctuations to 
generate noise can be reconciled, is by assuming phase 
cancellations from different regions of the jet. 

Plausible reason for the directional variation of 
correlation coefficients: It is known that the sound 
pressure fluctuations measured by a far field microphone 
is due to a sum of radiation from a large number of 
turbulent eddies distributed in the entire jet plume. 
Therefore, at first it was expected that individual point 
correlations, <p;p /> and <puu;p'>, would be very small. 
The experimental data, however, shows significantly high 
values at shallow angles and for the supersonic jet 
velocities. It is conjectured that present experiment picks 
up contribution from the coherent eddies and is unable to 
determine contribution from fine scale, small eddies. By 
definition, coherent eddies have long spatial coherence. 
The radiation from various spatial locations of such 
eddies are expected to be phase related (figure 18); 
therefore, a single point correlation reflects contribution 
from the entire eddy. Frequency analysis of cross- 
correlation data demonstrates that the measured 
correlations are from low Strouhal frequency 
fluctuations. Such fluctuations are known to be due to 
large organized structures which radiates primarily in the 
downstream direction close to the jet axis (Morris & 

Tam, 1979, Tam & Burton, 1984). Theoretical 
calculations (Michalke 1977, 1983) show that the higher 
the spatial correlation of turbulence, the narrower is the 
radiation angle. The quick drop in correlation with an 
increase of microphone polar angle, seen in figure 12, 
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further supports that the large spatially coherent 
structures have mostly contributed towards the measured 
correlations. Panda and Seasholtz (2002) have presented 
schlieren photographs of the same jet studied in the 
present paper. These photos demonstrate the presence of 
Mach wave emission (due to supersonic convective speed 
of some eddies) for the supersonic cases. The shock 
waves attached to the supersonically convected eddies 
makes an almost one to one connection between the near 
field turbulence and the far field sound pressure 
fluctuations. This is the reason for significantly higher 
correlation in supersonic plumes. The radiation fields of 
small turbulent eddies, however, is expected to be more 
omni-directional. In addition, they do not have significant 
coherence. Therefore, contribution from the small scale 
eddies towards single point <p;p'> & <puu;p'> 
correlations are expected to be extremely small, below 
experimental noise floor. This perhaps explains inability 
of the present scheme to identify noise sources for 90° 
radiation. 

IV. SUMMARY AND CONCLUSIONS: 

A recent advancement in molecular Rayleigh scattering 
technique, to simultaneously measure velocity and 
density fluctuations in high-speed flows, has been utilized 
to identify sound sources in supersonic and high subsonic 
unheated free jets. The particle- free, non-invasive 
technique involves passing a narrow CW laser beam 
through jet plumes and the collection of molecular 
scattered light from a point on the beam. A part of the 
collected light is directly measured to detect fluctuating 
intensity, which is related to air density fluctuations via 
calibration constants. The rest of the collected light was 
passed through a Fabry-Perot interferometer to detect 
Doppler shift associated with one component of jet 
velocity. The Fabry-Perot was operated in an imaging 
mode and the fringe formed at the image plane was split 
into two parts. The ratio of light intensity was related to 
velocity using another set of calibration constants. 

Two separate arrangements of collection optics 
were used to measure either p, u or p,v simultaneously. 
Time histories measured from various points in the plume 
were Fourier transformed to obtain spectra. A comparison 
shows that p and u spectra are similar in shape while v 
spectrum is different, especially in the centerline of the 
jet. The spectral data suffer from a large bias error due to 
electronic shot noise; yet the correlation data between the 
turbulent fluctuations and far field sound pressure 
fluctuations are relatively error-free. 

The noise emitted by the jets was measured by 
microphones placed in the far field and at polar angles 
from 30° to 90° to the downstream direction. The 
microphone signals p (t) were correlated separately with 
p(t), u(t), v(t), puu(t) and pvv(t) signals measured from 
various points in the flow. The non- intrusiveness of the 


laser-based technique avoids the probe-interference 
effects that have plagued previous attempts of source 
identification via causality method. Some significant 
observations from this study are the following: 

(a) Out of all flow parameters, the puu 
fluctuations are found to provide the highest correlation 
with the far field noise. This is closely followed by 
density p fluctuations. Indeed the <puu;p > and <p;p> 
correlations are strikingly similar in all respects: absolute 
magnitudes of normalized correlations, frequency 
distribution, and dependence on the microphone and 
probe volume locations. The v and pvv fluctuations are 
found to be poorly correlated with p . 

(b) The variation of the correlation coefficients 
with microphone polar angles, and frequency dependence 
of the cross-spectral data showed that the large scale 
organized structures mostly contributed towards the 
correlation coefficients. Cross-spectral analysis of puu 
and p fluctuations shows that the highest cross-spectral 
density is measured from Strouhal frequency, St ~ 0.2. In 
general, turbulent fluctuations in the frequency range 
0<St<l were correlated with far field noise. This Strouhal 
frequency range is typical of large organized structures. 

A study of correlation coefficients with microphone angle 
show that the polar angle closest to the jet axis (30°) 
provide the highest correlation; the magnitude falls 
sharply till 60°, beyond which data were mostly below 
experimental noise floor. Previous analysis (for example 
Michalke, 1977) has shown that the larger the spatial 
coherence of turbulence, the narrower the radiation angle. 
The long spatial coherence is characteristic of organized 
vortical motion. 

(c) The spatial dependence of the correlation 
coefficients was studied by moving the laser probe 
volume from point to point in the flow while keeping the 
microphone fixed (mostly at 30°). The strongest sound- 
producing region is found to lie along centerline and 
beyond the end of the potential core. For example, in 
Mach 1.4 jet the highest correlation was measured from 
centerline and x/D = 10, the potential core ended at x/D ~ 
7. When the laser probe was moved radially, the <pvv;p’> 
correlation was found to reduce monotonically and fall 
below the noise floor beyond the edge of the turbulent 
flow. Correlation from the lip shear layer is found to be 
Mach number dependent. Significant correlations were 
measured all along the lip shear layer of supersonic 
plumes, while the subsonic plumes did not show any such 
correlation. It is believed that the difference is due to the 
inception of the “Mach wave” emission process; where 
some eddies attain convective velocity higher than 
ambient. The far-field propagation of the shock waves, 
attached to the eddies, produce the high correlation in 
supersonic jets. In general, the higher normalized 
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correlation values measured in supersonic flow is also 
due to this Mach wave emission process. 

(d) To help various jet noise modeling efforts, 
Reynolds decomposition was performed on the full puu 
term, and the individual fluctuating terms were separately 
correlated with the far field noise. It was found that the 

first order fluctuations (2p uu'; p") and (p'uu; p") 
provided the largest contribution towards the full 
correlation <puu; p >. Terms which are second and third 
order in fluctuations contribute significantly less. This 
demonstrated, among other things, the limitation of 
traditional approximation, puu ~ puu, even in unheated 
subsonic jets. 

(e) According to the causality principle, second 
time derivatives of <puu; p > and <p; p > cross- 
correlations are related to the <p / ;p / > auto-correlation of 
far field sound pressure fluctuations. To avoid taking 
derivatives of experimental data such calculations were 
performed in frequency domain where cross-spectra are 
multiplied by freqency 2 . The principle dictates that a 
volume integral of all correlations over the plume leads to 
the far field sound pressure auto spectrum. In this work, a 
limited effort to determine effectiveness of a unit volume 
around a single measurement point to produce far field 
noise was calculated. It was found that the correlation 
values from the strongest noise producing regions are so 
high that a single point correlation is capable of 
producing more noise than the measured far field spectra. 
This has lead to a conjecture that phase cancellation from 
sound waves appearing from different parts of large 
coherent structures ultimately reduces the radiation 
effectiveness of a single source point. 
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Table I. Nominal Operating conditions: 


Sp. heat ratio, y = 1.4; Total temp., T 0 

= 300°K; Ambient density, p a 

= 1.16 Kg/m 3 , Ambient sound speed, a = 347 m/s 

Nozzle type 

Convergent operated at 
M = 0.95 

M 1.4 CD 

M 1.8 CD 

Minimum shock operation at M = 


1.395 

1.795 

Reynolds number Re D 

0.66xl0 6 

1.16xl0 6 

1.88xl0 6 

Jet Velocity Uj (m/s) 

316 

411 

486 

Estimated eddy convection speed U c 
(m/s) = 0.6Uj - 0.89Uj 

190-282 

247-366 

292-433 

Jet density pj (kg/m 3 ) 

1.36 

1.6 

1.89 

Frequency (kHz) for St = 1 

12.4 

16.2 

19.1 


Table II. Correlation coefficients between far field sound pressure fluctuations and various Reynolds decomposed terms 
of puu. 


Mach 

no. 

Prob 
e x/D 

Prob 
e r/D 

Mi 

c 

ang 

(puu;p') 

^puu'jp') 

(p'uujp 7 ) 

^upVjp^ 

(puV;p^ 

(pVu'jp') 

(puu) rms p m 

(pUU ) rms P rms 

(puuL.pL 

(puuL.pL 

(puuL.pL 

(p uu L s p ms 

1.8 

12 

0 

O 

O 

CO 

0.2207 

0.1635 

0.0636 

-0.00214 

-0.0025 

0.003 

1.8 

6 

0.45 

o 

O 

CO 

0.0572 

0.0406 

0.0201 

-0.0022 

-0.0023 

0.0011 

1.4 

10 

0 

o 

O 

CO 

0.1919 

0.1344 

0.0624 

-0.0021 

-0.0041 

0.0014 

0.95 

10 

0 

o 

O 

CO 

0.0682 

0.0592 

0.01 

0.0001 

-0.0019 

0.0007 

0.8 

8 

0 

o 

O 

CO 

0.022 

0.0198 

0.0037 

-0.0003 

-0.0012 

0.0003 
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Fig 1. Schematic of coordinate system and 
microphone locations. 



Fig 2. Principle of Rayleigh scattering technique. 


Images of the fiber face 




Fig 3. Fringes formed after passing through Fabry- 
Perot Interferometer by (a) incident laser light (b) 
Rayleigh scattered light. 



(a) (b) 

Fig 4. Principle of instantaneous velocity measurement by 
Rayleigh scattering technique. 


Microphone 



^microphone 


\ 

\ 



Fig 5. (a) Schematic of optical arrangement around jet facility to measure radial v component of velocity, (b) Scattering 
diagram for v measurement, (c) Top view of optical arrangement to measure axial u component, (d) Side view for u 
component measurement. 


NASA/TM— 2004-212385 


17 




(C) 



(d) 


Fig 5. (continued) (a) Schematic of optical arrangement around jet facility to measure radial v component of velocity, 
(b) Scattering diagram for v measurement, (c) Top view of optical arrangement to measure axial u component, (d) Side 
view for u component measurement. 


Auto-alignment 



Fig 6. Schematic of optical setup to analyze collected light. LI to L5 are lenses; BS1 is beam-splitter and FP1 is the 
Fabry-Perot interferometer. 
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Fig 7. Density and velocity calibration using 3 photo-multiplier tubes. 



probe at indicated locations of Mach 0.6 plume. 
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Fig 9. Power spectral density of (a) air density, (b) axial velocity, (c) radial velocity, (d) density*(axial velocity) 2 
fluctuations; (e) density-axial velocity cross-spectrum, and (f) phase of cross-spectrum in Mach 1.4 jet at centerline & 
x/D = 10 
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Fig 10. Normalized cross-correlation between sound pressure fluctuations p / and (a) air density p, (b) puu, (c) axial 
velocity u, (d) radial velocity v & (e)pvv fluctuations measured at x/D=10 and centerline of M=1.4 jet. The microphone 
was kept at 50D and 30° to jet axis unless otherwise mentioned. 
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Fig. 11. Normalized cross-correlation between sound pressure fluctuations and (a) puu fluctuations (b)-(f) various 
Reynolds decomposed terms of puu measured at x/D=10 and centerline of M=1.4 jet. The microphone was kept at 50D 
and 30° to jet axis. 
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Fig. 12. Angular dependence of peak correlation between microphone pressures and (a) p, (b) puu fluctuations. The 
laser probe was kept fixed at centerline & x/D=8 (M=0.95 case) or 10 (M=1.4), or 12 (M=1.8) while the microphone 
was moved at various angular locations on a 5 ID arc. 



Fig. 13. Axial dependence of peak correlation between microphone pressures and (a) p, (b) puu fluctuations 
measured along centerline of indicated Mach number plumes. Microphone was kept fixed at 50D and 30° to 
jet axis. 
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Fig. 14. Axial dependence of peak correlation between microphone pressures and (a) p, (b) puu fluctuations 
measured along shear layer(r/D=0.45) of indicated Mach number plumes. Microphone was kept fixed at 50D 
and 30° to jet axis. 



Fig. 15. Radial dependence of peak correlation between microphone pressures and (a) p, (b) puu fluctuations, measured 
at x/D=10 and various radial positions, in Mach 1.4 plume. Microphone was kept fixed at 50D and 30° to jet axis. 
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Fig. 16. Normalized cross-spectrum of sound pressure fluctuations p and (al to el) air density fluctuations, 
(a2 to e2) puu fluctuations from indicated probe locations in various Mach number plumes. Microphone was 
kept at 30° and 50D away. 
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Fig. 17. Relative contribution to the 30°-microphone auto-spectrum from pvv-p / and p-p / cross-spectra measured in 
Mach 1.4 plume; probe location x/D=10 & centerline, (a) p-p' cross-spectra, (b) pw-p' cross-spectra. Calculated 
contribution from (c) dipole (equation 14), and (d) quadrupole (equation 13) terms. 



Fig. 18. Schematic of large and small scale turbulent structures and noise radiation. 
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